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Abstract 

A Direct Numerical Simulation of a coolant jet injected normally into a periodic crossflow 
is presented. The physical situation simulated represents a periodic module in a coolant hole 
anay with a heated crossflow. A collocated finite difference scheme is used which is fifth-order 
accurate spatially and second-order accurate temporally. The scheme is based on a fractional step 
approach and requires the solution of a pressure-Poisson equation. The simulations are obtained 
for a blowing ratio of 0.25 and a channel Reynolds number of 5600. The simulations reveal the 
dynamics of several large scale structures including the Counter-rotating Vortex Pair (CVP), the 
horse-shoe vortex, the shear layer vortex, the wall vortex and the wake vortex. The origins and the 
interactions of these vortical structures are identified and explored. Also presented are the 
turbulence statistics and how they relate to the flow structures. 

Introduction 

In recent years, Direct Numerical Simulation (DNS) has become a powerful tool in 
understanding flow physics and the interaction between various dynamical scales in the flow. 
Numerical Simulation of channel flows, due to its geometrical simplicity, has been studied most 
commonly in the DNS literature. Orszag et al. (1980) studied transition from laminar to turbulent 
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flow using a 3-step time splitting method. Fourier series expansions in the spanwise and 
streamwise directions and Chebyshev polynomial expansion in the direction normal to the wslls. 
Orszag et al. (1983) later used the same algorithm, but with an explicit Adams-Bashforth 
integration scheme, to perform a DNS of turbulent channel flow. A DNS for a curved channel 
flow was carried out by Moser et al. (1984) at Re. = 1 80 again using a spectral method. Kim et al. 

(1987) used a Fourier Chebyshev pseudospectral method to study fully turbulent channel flow at 
the same Reynolds number and presented a large number of turbulence statistics. A finite 
difference DNS for a turbulent channel flow was presented by Rai et al. (1989) who showed the 
need for higher order accurate upwind schemes to yield accurate simulations. 

These initial successes with the DNS method spurred the investigation of slightly more 
complex situations. Gavrilakis (1992) and Huser and Biringen (1993) both reported DNS for flow 
in a square duct geometry. Lyons et al. (1991) studied fully developed passive heat transfer in a 
direct simulation based on spectral method. Kasagi et al.(1992) and Kim and Moin (1989) 
obtained the solution to the energy equation for channel flows using a DNS procedure. Sumitani 
and Kasagi (1995) have reported a DNS in a channel flow with uniform suction and injection. 
More recently, Pointel et al (1996) have presented a DNS of channel flow with heat transfer 
using a collocated fifth-order accurate finite difference procedure. The calculations were done for 
a Re,= 180 on a 128x64x64 grid and were found to be in good agreement with the DNS data of 
Kim et al. (1987) and the measurements of Kreplin and Eckelmann (1979). 

The present paper aims to extend the work of Pointel et al. (1996) for a more complex 
flow situation, and intends to perform a DNS of a periodic cross-flow in a channel with a vertical 
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slot-injection of coolant through the bottom surface. Examples of slot injection into a crossflow 
include the film cooling of gas turbine blades and combustor walls, flow control of aircraft and 
missiles through blowing, and effluent discharges into air or water bodies. The application of 
primary interest here is the film cooling of gas turbine blades. 

Jets in crossflow, as encountered in film cooling configurations, represent a rather 
complex flow situation. Several investigators have studied the features of the flow field 
associated with a single jet in a crossflow. The counter rotating vortex pair (CVP) is the most 
dominant structure observed, and persists far downstream of the jet injection. The CVP is a 
consequence of the reorientation of the jet-hole vorticity by the crossflow (Foss (1980), 
Andreopoulos (1985), Kelso et al. (1996)). Shear layer or ring vortices have also been observed in 
the near field of the jet. These vortices are driven by the Kelvin-Helmoltz instability. Horseshoe 
vortices are produced as a consequence of the jet obstructing the crossflow producing an adverse 
pressure gradient just upstream of the jet. These vortices have spanwise vorticity when they 
form, and as they travel downstream they are deflected by the crossflow and acquire stream wise 
vorticity rotating in a direction opposite to the CVP. Wake vortices which exist between the wall 
and the jet have also been observed (Fric and Roshko, 1994). 

The majority of the studies aimed at numerically examining the details of the flow 
structures have approached the problem by solving the Reynolds Averaged Navier Stokes 
(RANS) equations. This time-averaged approach requires the use of turbulence models and makes 
predictions of the dynamic nature of vortices unfeasible. RANS calculations have been reported 
by Patankar et al. (1977), Sykes et al. (1986), Kim and Benson (1992) and Demuren (1993). 
Recently Large Eddy Simulations (LES) that resolve the dynamics of the large scales only and 
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model the small scales have been presented by Jones (1996) and Yuan and Street (1996). 

The geometrical configuration of interest considered in this paper is shown in Figure 1 . 
The primary flow (called the crossflow) in the channel is in the x-direction. The coolant jet is 
injected vertically upwards from a rectangular slot on the bottom surface. Periodic boundary 
conditions are imposed on the crossflow in the x-direction. The periodic flow in the x-direction 
makes the computational problem not only simpler to handle but it also allows us to mimic the 
effects of upstream film coolant holes. In a film cooling configuration with a series of parallel 
holes in the streamwise direction, the downstream coolant hole will see the downwash from the 
upstream coolant hole, and this effect is incorporated in the present configuration with the 
streamwise periodic boundary condition. Single coolant hole studies with Dirichlet inflow and 
outflow conditions can not simulate the upstream hole effects. The correct alternative is to 
include all rows of holes in the calculation domain, but performing a DNS calculation for such a 
configuration is considerably more difficult than the simpler single-hole situation considered in 
this paper. However, for a single hole domain, the imposition of periodic boundaries in the 
crossflow direction is inappropriate, unless a suction boundary condition is incorporated, since 
with injection alone there is a net addition of mass across the lower surface. Therefore, in order 
to impose periodic conditions in the primary flow direction, a slot-suction boundary condition is 
imposed on the top wall. The normal suction velocity and temperature at the top wall is made 
exactly equal to the normal injection velocity and temperature at the bottom wall so that the 
mass, momentum and energy injected in exactly balances that sucked out. The blowing (or 
suction) ratio (injection velocity/crossflow velocity) is taken to be sufficently small (0.25) so that 
the injection and suction effects are primarily confined to a small distance normal to the surface. 
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In the span wise direction, a row of film coolant holes are assumed to be present. 
Therefore, in the computational domain of this study, consisting of a single hole, periodic 
spanwise boundary conditions are imposed. 

Numerical Procedure 

Assuming constant density, and using the channel half-width 5 as the length scale, and the 

friction velocity u x for the fully-developed channel flow (with no injection or suction) as the 

velocity scale, the non-dimensional continuity equation and the Navier Stokes equation for each 
velocity component u, (i=l, 2 and 3 represent the u, v, w velcoity components) can be written as: 

V.u=0 0) 

— + (u.V)u, = -Vp + j^-V 2 u i (2) 

dt T 

In the above equation, p represents the pressure and Re t represents the Reynolds number based 

on the fully-developed channel flow friction velocity. For the temperature, iso-flux, iso-thermal 
conditions proposed by Kasagi et al. (1992) have been adopted. According to their investigations, 
for iso-flux conditions, the variations in the wall temperature for air flow can be considered to be 
negligible. Hence we can assume the wall to be locally isothermal. The isoflux condition implies 
that the ensemble averaged wall temperature (over z and t), increases linearly with x-direction. 
Hence the bulk mean temperature averaged over z and t also increases linearly in x-direction and 
can be obtained from an energy balance. We also assume that the injection and suction jets are at 
the same temperature but lower than that of the wall. Hence the overall energy balance is 
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preserved. The non-dimensional energy equation can then be expressed as: 


— + (mV)0 
di 


Re T Pr 


v 2 e 



( 3 ) 


where 


1 X 


(4) 


and T x is the friction temperature given by q^(pc p u t ), and q w and C^w) are the ensemble averaged 
heat flux and wall temperature (over z and t) respectively. 

Results are obtained for a Re t =180. This corresponds to a channel Reynolds number Re m 
based on the mean velocity and the channel height of 5600. As noted earlier, DNS results for 
fully developed channel flow at this Re t value were reported by Kim et al. (1987) and Pointel et 
al. (1996). Since in the configuration considered in this study, the net mass addition is zero, 
Re t =180 represents an average value for the channel. The injection and suction jet-holes were 

placed in the center of the lower and upper walls and were rectangular in shape with the long 
dimension of the hole (d) along the crossflow direction (x) and the short dimension of the hole 
(w) in the spanwise direction (z). The corresponding non-dimensional hole sizes d/8 and w/5 
were 0.5 and 0.2 respectively. The actual non-dimensional lengths of the computational domain 
were chosen to be 4n x 2 x 0.4*. These correspond with the 4n x 2 x 2n computational domain 
dimensions of Pointel et al. (1996), except in the spanwise direction where coolant hole spacings 
of 2.5d (=1 .25 = OAn in non-dimensional units) has been assumed. 
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As noted earlier, the flow and temperature variables are assumed to be periodic in the x 
and z directions. Dirichlet. zero velocity and temperature, boundary conditions were imposed at 
the channel walls while a zero-value Neumann condition was imposed on the normal pressure 
gradient at the wall. At the injection and suction slots, the normal velocity is set equal to 0.25 of 
the mean crossflow velocity corresponding to the fully developed channel flow' at Re =180, and 
the temperature is set to a value lower than the wall temperature, and is assigned a non- 
dimensional value of 0.5. 

The Fractional Step approach originally developed by Chorin (1967) was used for 
integrating the equations. In the first step, the momentum equation without the pressure gradient 
contribution is solved using a second order explicit Adams-Bashforth time-integration for the 
convection and streamise and spanwise diffusion terms. For the normal-diffusion term, an 
implicit Crank-Nicolson scheme is used. This approach uses an intermediate velocity (denoted 
by a ' ), and for a typical velocity component, say u, has the following form: 

llill = i(C" + Dr ) - + A”"' ) + A" + A )• <4) 

Al 2 2 l 

Here C represent the convection terms, D e represents the streamwise and spanwise diffusion 
terms and Dj represents the normal diffusive term. 


The velocity field is computed in the second step: 



where the pseudo pressure <t> is related to the pressure by: 


We can apply the divergence operator on equation (5) and use continuity of velocity field to 

obtain the pressure Poisson equation: 

-LvJ7 = V.(V^'). (7) 

A/ 

Once the pseudo pressure is known the velocity field u n+1 can be calculated from equation (5). 

The scalar field for temperature is then obtained by integration of the energy equation. We 
use a similar approach to the momentum equation to get the fractional step equations from which 
we can directly obtain the temperature field. 

The spatial discretization was done using the non-conservative formulation of the 
convective and diffusive terms on a collocated grid with velocity and pressures at the center of 
the control volume. The convective terms were treated by a fifth order upwind biased 
approximation as used by Rai et al. (1989). The diffusive terms are treated by a sixth order 
accurate seven point central difference stencil. For the boundary terms we used one sided 
differences to maintain the accuracy of the scheme. The pressure Poisson equation is discretized 
using the 4-2 consistent formulation as recommended by Najjar et al. (1994). This involves 
second order discretization of the divergence operator and fourth order discretization of the 
pressure gradient operator. In the continuity equation, which is used for updating the velocity 
field at the cell faces, the gradient operator is represented by a fourth order finite difference 

formulation. 
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The implementation of the numerical procedure, described below, is done on a MasPar 
MP-2. This is a massively parallel computer and provides a suitable platform for computer 
intensive DNS calculations. The MasPar is a SIMD machine, which has a 128 x 64 parallel 
processors (PE) array. We have used a grid of 128x64x64 for the calculation, which is uniform in 
the x and z direction, and non-uniform in the y direction. Each (I, J) node is mapped on to a PE. 
Each PE has a local memory of 64K and provides the storage in the third (K) dimension. 

The time integration was first done for the velocity field till the flow reached a fully 
developed state. Then the thermal field was allowed to reach the fully developed state before the 
calculation of mean quantities was started. Typically 60,000 time steps were required to get to 
this stage. The calculation of the mean quantities was done over 20,000 non-dimensional time 
steps. This was followed by the calculation for the turbulent statistics and this was done over 
8000 steps. During the statistics calculation we also calculated the instantaneous vorticity 
profiles for various streamwise, crossstream and spanwise planes at various time steps m order 
to present the time evolving characteristics for the flow field. 

Results and Discussion 

In making an assessment of the adequacy of the grid size used in the present study where 
the mean Re t is 180, we make comparisons with the grid resolution of reported channel flow 
studies at a Re t of 180. Pointel et al. (1996) used a 128x64x64 grid for simulating fully developed 
channel flow in a 4n x 2 x 2rt domain at a Re t of 1 80. Results obtained are shown in Figure 2 and 
compared with the DNS predictions of Kim et al. (1987) and the experimental data of Kreplin 
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and Eckelmann (1979). As can be seen, the comparison is quite good. The same 128x64x64 grid 
was employed in the present problem which has a smaller computational domain of 4n x 2 x 

0.47t, and therefore has an improved resolution. The resolution in wall-coordinates are noted to 
be: Ax* = 17.67, Ay* = 0.018 - 23.76, Az* = 3.59. In comparison, Kasagi et al. (1992) reported 
corresponding grid spacings of 18.4, 0.08-4.9 and 7.4 in a 5 ti x 2 x 2n computational domain. 
Similar grid spacings of 17.7, 0.05-4.4, and 5.9 were also used by Kim and Moin (1989) and Kim 
et al. (1987) in a 4 tc x 2 x 2n computational domain at the same Reynolds number. Moser and 

Moin (1984) used 128x65x128 points in a similar computational domain and comparable 
Reynolds number, and reported that this resolution was adequate to resolve all the essential 
turbulent scales. Lu and Hetsroni (1995), in their direct simulation of channel flow with heat 
transfer at a comparable Reynolds number, have used Ax*, Ay*, and Az* values of 32, 0.4-20, 
and 8. Based on these comparisons, the grid size used in this study should be deemed to be 
adequate. 

In presenting the results below, the origin of the coordinate system is fixed at the jet-hole 
center. The x and y dimensions are scaled with the long dimension of the jet hole (d) oriented in 
the direction of the crossflow, while the z dimension is scaled with the spanwise dimension of 
the jet hole (w). 

Characteristics of the Flowfield 

Figures 3(a) and (b) present the u and v velocity contours in the spanwise midplane (z/w 
= 0), and correspond to the vector plots shown in figures 5(a) and 5(b). For clarity, several of the 
plots are stretched in the vertical or spanwise directions. Figure 3(a) and 5(a) clearly show the 
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almost immediate bending of the jet due to the strong crossflow. The jet penetration is confined 
to a y/5 of about 0.5 or a y/d of about 1 which corresponds to the long-dimension of the jet-hole 

(hereinafter referred to as the jet diameter). Figure 5a clearly shows that at the channel mid-plane 
(y/d=2), the crossflow velocity is relatively unaffected by the injection or suction slots. This 
observation supports the earlier premise that the behavior of the film-coolant injection jet can be 
studied independent of the suction jet, and that the observations made should hold for a film 
cooling configuration with a series of coolant holes in the crossflow (or primary flow) direction. 

A weak horseshoe vortex system is observed in Fig. 3(a) and in the near-wall detail of the 
velocity vectors shown in Fig. 5(b). The vortex system extends over one hole diameter upstream 
of the jet and appears to consist of two clockwise eddies rotating side by side. However, the 
horse-shoe vortex system does not appear to be deflected around the jet as it does at higher 
blowing ratios (see Fric and Roshko, 1994; Kelso et al., 1996; Muldoon and Acharya, 1998). No 
footprint of the vortex is seen in the x-vorticity contour along a horizontal plane just above the 
lower surface (Fig. 3c). 

In examining the instantaneous v- velocity contours (Fig. 3b), it can be seen that roughly 5 
hole diameters downstream of the injection location (x/d=5), alternate pockets of positive and 
negative v-contours are present (note the high level of stretching in the y-direction). These are 
analogous to those observed by Kim et al. (1987), and represents ejection and sweep events 
associated with a developing boundary layer. In the near vicinity of the injection hole, 1 < x/d < 
5, the behavior is different with primarily positive near wall velocities. 

The x-vorticity contours shown in Figure 3c indicate the existence of a dominant wall- 
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vortex structure close to the lower surface with vorticity opposite to the CVP. Thus downstream 
of the injection hole, positive vorticity is seen close to the wall in the right half (viewing the flow 
from a downstream point) and negative vorticity is noted in the left half. Further away from the 
wall, at a higher y/d value, the CVP is encountered and the sign of the vorticity reverses. This can 
be seen in the cross-sectional (y-z) plots of x-vorticity presented later. Figure 3c also shows the 
generation of negative and positive x-vorticity on the right and left edges of the jet hole; these 
represent the source of vorticity for the CVP. It is however important to note, that the vorticity 
for the CVP does not originate from any in-hole boundary layer vorticity (as has often been 
reported in the literature), since at the jet exit plane a uniform v-velocity was specified. Rather, 
the x-vorticity along the jet-hole edges is generated by the interaction of the injected velocity and 
the crossflow. 

Figure 4 shows the velocity vectors in different cross-stream planes at x/d=0.5, 2, 5 and 
10. At x/d = 0.5 we can see the jet penetration to be more than y/d =1.8. However beyond a y/d 
of 1, the velocity magnitudes in the y-z plane are considerably smaller than the velocities in the x- 
y plane, and therefore as seen in Fig. 5a, the flowfield is relatively unaffected by the jet beyond a 
y/d of 1. We can see the evolution of the CVP around z/w = ± 0.8 for y/d < 0.2. There is 
considerable flow across the spanwise periodic boundaries indicating a significant level of lateral 
transport, and interaction between spanwise adjacent coolant jets. Both inflow and outflow are 
noted across the periodic boundaries. We also note the absence of any symmetry across the 
z/w=0 boundary. In fact by x/d=10, the CVP has lost its bound- vortex structure and no 
semblance of symmetry exists across z/w=0. In the near field (x/d<5), each half of the CVP is 
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clearly associated with a weak counterrotating eddy system above it. 

Figure 5c shows that there is a significant amount of entrainment of the crossflow into the 
wake of the jet. This clearly appears to persist up to about x/5 of 9 at which point the wake 
appears to meander away from the jet-hole centerline. A close inspection also reveals flow 
reversals of the entrained crossflow near the jet-hole centerline. These reversals are most clearly 
seen between x/d of 3 and 5 in Fig. 5c, and represent the onset of the wall-vortex system which 
will be clearly evident in the x-vorticity contours in a y-z plane presented later. 

Figure 6 shows the profiles of the mean velocity and temperature, non-dimensionalized 
by the corresponding average friction-velocity u T (which is the same as the friction velocity for 
fully developed turbulent channel flow at the same average Reynolds number Re t =180 of this 
study) and the average friction-temperaure T t . The profiles are presented for several x/d locations 
and two z/w planes, z/w=0 corresponding to the jet-center plane, and z/w=-1.9 representing a 
location between the jet-hole and the spanwise periodic boundary. The near-wall behavior 
(y + =y.Re/6 < 12), appears to conform somewhat to the law-of-the-wall behavior with a 

maximum deviation of less than 30%. However, in the turbulent region (y + >12), log-law behavior 

is not observed even off the jet-center plane and upstream of the jet (z/w=-1.9, x/d=-2) due to the 
effect of vortical flow structures. At z/w=-1.9, the near- wall behavior appears to be only mildly 
altered by the spanwise cross-stream flow downstream of the jet In fact. Fig. 5c indicates that 
the spanwise cross-stream flow close to the wall at z/w-- 1.9 is relatively weak. The maximum 
deviation from the x/d=-2 profile appears to occur in the vicinity of x/d=2, beyond which flow 
recovery is initiated. The profiles in the jet-center plane (Fig. 6b and 6d) indicate significant 
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distortion of the near-wall behavior by the jet This distortion begins at the upstream end of the 
jet (x'd=-0.5) with a significant reduction in the wall-gradients. Flow recovery begins downstream 
of the jet. but the recovery in the inner-region (y*<12) is substantially faster than in the outer- 

region (y*>12), leading to a change in the sign of the slope in the velocity profile in the buffer 
region and immediately above it (7<y*<35 or 0.08<y/d<0.39). Such changes in the slope of the 

velocity profile can lead to the development of flow instabilities and large-scale vortical 
structures. This y/d region can be seen to correlate with the region associated with the growth of 
large scale structures in the vorticity contours presented later. 

Figure 7 shows the time history of spanwise vorticity at z/w = 0. The time-series is 
presented at time instances t|, t 2 ,...ts, separated by At=0.02 where At is the non-dimensional time 

interval. The figures indicate the predominant negative to z on the windward side which is 
indicative of the crossflow bending due to the jet. Patches of positive vorticity can be seen on the 
leeward side of the jet. The vortices on the windward and leeward side of the jet are driven by 
Kelvin-Helmholtz instabilities, and are generally referred to as shear layer or spanwise-jet- 
vortices. These vortices, particularly those from the windward side, appear to be peeled by the 
crossflow and convected downstream. In fact the near-wall boundary-layer region downstream of 
the jet is dominated by the negative vorticity partly originating from the windward side of the jet, 
and partly associated with a developing boundary layer. The dynamical nature of the vortices 
should be noted in the time-series presented. Note the vortices on the windward side at time t 3 
between y/d of 0.5 and 1. At the next time-instance, a necking process is initiated, and 
subsequently two vortices form and convect downstream. A third vortex is transported to the 
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near wall region downstream of the jet. Upstream of the jet the horseshoe vortex system 
consisting of two eddies can be observ ed, but it can be seen to be rather small and localized. 

Figure 8 shows the instantaneous x-vorticity at x/d=-0.5 (upstream edge of the injection 
hole) and at x/d=0.5 (downstream edge of the injection hole). At x/d=-0.5, the begimngs of the 
CVP can clearly be seen. This was also observed in Fig. 3c, and implies that the x-vorticity of the 
same sign as the CVP is originated on the windward side through the interaction between the jet 
and the crossflow. Both at x/d=-0.5 and 0.5, the CVP appears to be straddled by an eddy rotating 
in the opposite direction. This is possibly the horse-shoe vortex that lifts up above the CVP at 
this low blowing ratio. Reported studies in the literature at higher blowing ratios indicate that the 
horse-shoe vortex is deflected in the spanwise direction by the jet and travels downstream along 
the wall parallel to the CVP. In the present case the signature of the horse-shoe vortex in the 
spanwise direction is not observed, and we believe that at the low blowing ratio considered here, 
since the jet is immediately bent by the crossflow, the horse-shoe vortex instead of being pushed 
in the spanwise direction is pushed up in the vertical direction to ride above the CVP. 

Figure 9(a) shows the time-series (presented 2At apart) of the x-vorticity at x/d=2. Three 

distinct structures can be observed: the CVP, the horse-shoe vortex system, and a wall vortex 
below the CVP. The orientation of the wall vortex and the horse-shoe vortex are in the same 
direction and opposite to the CVP. The wall vortex is due to the entrainment of the crossflow 
into the wake of the jet. The entrained crossflow encounters high pressures along the jet-hole 
centerplane which drives the flow along the wall away from the centerplane, thus establishing a 
wall-vortex. The horse-shoe vortex is now shifted outwards toward the periodic boundaries, but 
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still remains lifted off the surface. In examining the time-series, one notes that the CVP undergoes 
significant distortion with time. There is also significant spanwise transport of vorticity across 
the periodic boundaries leading to fairly effective spanwise mixing. 

At x/d=5 (Fig. 9(b)), a noticeable interaction between the wall-vortex and the CVP is seen. 
The wall-vortex can be seen to be entrained upwards towards the CVP, and in this process is 
reoriented into y-vorticity representing wake vortices (seen in Fig. 10 and 1 1). This reorientation 
of the wall-vortex system into a w ake vortex system has also been reported by Fric and Roshko 
(1994) and Kelso et al. (1996). 

Further downstream at x/d = 10 (not shown), we cannot clearly identify the CVP as 
being an organized structure spanning both sides of the jet-hole vertical midplane. There is 
significant vortical activity though, but of a much more complex nature, and indicates the 
persistence of the wake. 

Figure 10 shows the instantaneous (Dy contours close to the lower surface at y/d = 0.002, 
while Fig. 11 shows the time history of cOy contours at y/d = 0.0239 respectively. In Figure 10, 
we note that the wake vortices strongly manifest themselves primarily beyond x/d=5 which is the 
approximate location where we noted upward entrainment and reorientation of the wall vortex 
system (Figure 9(b)). This observation provides further evidence of the correlation between the 
wall-vortex and the wake-vortex system. The wake vortices are in the form of a series of co- 
rotating or counter-rotating adajacent eddies that have features quite distinct from those behind a 
rigid cylinder. At earlier x/d locations, upstream of x/d=5, the wall vortex had remained fairly 
stable and bound to the near wall-region, and no significant presence of a series of adjacent co- 
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and counter- rotating wake vortices was noted. Instead, y-vorticity in the near wake (x/d<5) was 
observed in the form of a near-contiguous elongated chain (see Figure 10 and 11) and stems form 
the interaction of the jet and the deflected crossflow which leads to significant values of the du/dz 
term (and therefore «,) along the spanwise edges of the jet. This vorticity is entrained into the 
near wake (x/d<5), and eventually breaks down further downstream, contributing also to the 
strong wake vortex system observed beyond x/d of 5. The time series behavior shown in Fig. 1 1 
also indicates that the wake structure oscillates around the jet-hole center-plane in the near field, 
but exhibits significant asymmetry in the far field. At time instances t h t 3 and t 5 , for example, the 
wake vortex system beyond x/d of 6 is skewed considerably towards positive z/w values. This 
departure from symmetry was also noted earlier in velocity vector plots of Fig. 5c. 

Figure 12 presents the instantaneous pressure contours at y/d=0.0239 at four time- 
instances. The adverse pressure gradient upstream of the jet can clearly be seen, and is 
responsible for the horse-shoe vortex system. In the jet-hole region, and downstream of it, the 
pressure distribution shows a necklace pattern. Note that in the region -0.5<x/d<0.5 there is a 
strong pressure gradient in the crossflow directed inwards in the spanwise direction toward the 
jet. This inward pressure gradient combined with the adverse pressure gradient the flow 
encounters at the jet interface leads to the generation of x-vorticity of the same-sign as the CVP 

along the spanwise boundaries of the jet-hole (z/w= ±0.5). This mechanism therefore appears to 

be responsible for the initial development of the CVP. 

Downstream of the jet-hole, the spanwise pressure gradient is directed inwards towards 

the jet centerline. This pressure gradient is responsible for the inward entrainment of the 
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crossflow that is the source of the wall vortex system. Note that downstream of about x/d-5, the 
coherence of the pressure contours is somewhat lost leading to the asymmetric wake pattern 
reported earlier. 

Time-averaged Statistics 

Figure 13 and 14 show the rms contours at x/d=0.5, 2, 5 and 10. At x/d=0.5, and w ms 
are comparable in magnitude with the latter attaining slightly larger values, but u^ is nearly twice 
the Vnnj and w^ values. The u^ and w^ plots correlate well with each other exhibiting maximas 
in regions of high unsteady (0* vorticity, that is in the outer regions of the CVP and the CVP- 
horseshoe vortex interface (see Fig. 8) that exhibit a significant level of unsteadiness, and near the 
spanwise periodic boundaries across which significant vorticity transport occurs (Fig. 8 shows 
only one instance in time, and shows high vorticity near the z/w=-3 boundary . At other instances 
in time high vorticity is seen at the z/w=3 boundary as shown in Fig. 9). The v^ and t m5 profiles 
correlate well with each other, but do not show a strong correlation with the u ms and w„ plots. 
In fact the peak values are obtained near y/d=l, and appears to correspond to the extent of the 
vertical jet penetration. The interface between the jet and the crossflow, is characterized by 
Tollmein-Schlichting instabilities and high spanwise vorticity <»* (see Fig. 7), and is therefore 
associated with high v^. Since the jet and the crossflow are at different temperatures, their 
interface is also associated with high W The aforedescribed correlation between the rms values 
and the coherent structures imply that the large scales are largely responsible for the production 
of turbulent stresses. At x/d=2, the behavior of the rms profiles is quite similar to that at x/d=0.5, 
except that the W profiles shows high values both in regions of high v^ (interface between the 


TT1 


18 


upper edge of the jet and the crossflow) and high u™ (interface between the CVP and the 
crossflow). Further downstream (x/d=5 and 10) the jet deflection is complete, and the CVP 
dominates growing in size and strength. All the rms-profiles now correlate with the high vorticity 
regions of the CVP, and as a consequence correlate well with each other. At x/d=10, we still 
observe Urms > v ms > w^. although the relative magnitudes are closer to each other than they were 
at x/d-0.5. Note that the time-averaged rms-profiles at x/d=10 show reasonable symmetry 
around z/w=0, while the instantaneous profiles shown earlier (Fig. 5c and 8-11) displayed 
significant asymmetry in the spanwise direction. 

The rms profiles in the z/w=0 plane are shown in Figure 15. All the peak values occur 
downstream of the jet injection (x/d=0), but the upwash effect is clearly evident. High rms values 
enter the computational domain from the upstream boundary, and the flow development 
downstream of the injection point is constrained by these high rms values. The u^ value peaks 
first in the x-direction (near x/d=5) and occurs closest to the wall (y/d=0.25), while the v^ peaks 
around x/d=6.5. These peak locations correspond with the location where the shear layer vortices 
saturate in strength (see Fig.7). The other rms values and W peak at an x/d in the vicinity of 
7. Note that in the vicinity of this x/d (in the range of 5-7), significant reorientation of the wall 
vortex system into a wake vortex system occurs, and therefore the wake- and wall-vortex 
structures structures apparently play an important role in the and W values. 

Figure 16 and 17 present the Reynolds stress components at x/d=0.5, 2, 5 and 10. As 
expected, uV is observed to be symmetric about z/w = 0, and u'w' and VW are observed to be 
antisymmetric about z/w = 0. The high uV values are associated with the jet-crossflow interface 
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shown in Fig. 7. where shear layer vortices primarily form. At x/d-10, distinct regions of upper 
and lower interfaces are no longer clearly ev ident and a large patch of high uV is noted in the jet- 

flow region. High values of uV are also associated with the vortical structures transported 
across the periodic boundaries. These structures lead to large values of 3u/dy and therefore uV. 
The u'w' and v'w' values are high in regions of large du/dz and dv/dz respectively. These regions 
are again associated with the jet-crossflow interfaces and the vortical structures crossing the 
periodic boundaries. Note that the clover-leaf shape and the sign of the v'w' contours near the 

wall for x/d > 5 indicate the existence of counter-rotating eddies near the wall, i.e., the CVP and 
the wall vortex, in the y-z plane. In comparing the magnitudes of the shear stresses, in the near 
field (x/d < 2), u'w' > uV > v'w' with u'w' being substantially higher presumably due to large 

BuJdz values near the edges of the jet. Further downstream, u'w' and u'v become more 
comparable in magitude. 

Figure 1 8 presents the Reynolds stresses in the z/w = 0 plane. The peak shear stresses, 
in all cases occur downstream of x/d=5. The u'v' contours correlate best with the v^ contours 
(Fig. 15) and therefore correlate with the shear layer vortices, and peak at the locations where 
these vortices saturate in strength. At x/d=5, for example, the peak uV value shown in Fig. 1 7 is 
2.85 and its y/d location correlates with the shear layer vortices, while secondary bimodal peaks 
(of the order of 1) occur at y/d=0.25 which correlates with the CVP and the horseshoe vortices. 
The v'w' contours appear to be consistent with the nature of the near-wall flow structure in the 

y-z plane. Close to the wall, the wall vortex produces negative v'w/ values. The wall vortex 
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system gains prominence, as seen earlier, around x/d-5 and therefore large negative v w values 
are obtained in the vicinity of this x/d value. Further downstream, the wall vortex is entrained as 
an upflow into the CVP. and as seen in Fig. 1 1, wake vortices are noted. The v'w' contours in 

this region show a more complex behavior and are positive. The u w contours also correlate well 
with the crossflow entrainment. They peak where the crossflow is a maximum (around x/d-5), 
and decay further downstream. 

Figure 19 and 20 present the turbulent heat flux correlations at x/d=0.5, 2, 5 and 10. The 
turbulent heat flux correlations appear to correlate reasonably well with the rms values and the 
behavior of the large scale vortical structures, thus indicating the important role these structures 
play in the turbulent heat transport. In comparing the magnitudes uY > wY > vY in the near 
field (x/d<5) indicating that the stream wise velocity fluctuations are most responsible for the heat 
transfer from the surface. Further downstream, at x/d=10, the vY values become comparable in 
magnitude to the uY values, while the wY decrease in magnitude and become less important. 
The heat flux distributions at z/d=0 in the x-y plane are shown in Fig. 21, and support the above 
observations. The maximum wY occurs near x/d=5 where there is significant spanwise 

entrainment of the crossflow. The peak in the uY and vY occur downstream of x/d=5 and at a 
y/d location of approximately 0.5. These peak values are therefore associated with the spanwise 
shear layer vortical structures (Fig.7). 

Figure 22 presents the turbulent viscosity factor v tii calculated as: 
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where all terms on the right hand side of the above equation are computed from the DNS results. 
The viscosity factors are presented in a x-z plane corresponding to a y/d=0.1085. While all three 
viscosity' factors ^t.12' ^t.13’ ^t. 2 3 were computed, only is presented. Note that the 
magnitudes of the individual viscosity factors ^t,i2» W,i3> W,23 at the same point are 
substantially different in magnitude indicating the anisotropy in the flowfield. The ^t,i2 clearly 
corresponds to the wake structure shown in Fig. 11, and indicates that the large-scales clearly 
have a role in the values of the turbulent viscosity. The footprint of the horseshoe can also be 
seen upstream of the jet hole. The ^t,i 3 * v tj 23 values show less coherence with respect to the 
large scale structures, and are characterized more by patches of high positive and negative 
viscosity factors. The positive and negative values of the viscosity factors are associated with the 
fact that both the turbulent shear stress and the mean strain rate can have both positive and 
negative values. Clearly, any use of an isotropic turbulent viscosity is going to be unsuccessful in 
capturing the true behavior. 

Concluding Remarks 

A Direct Numerical Simulation of a film coolant jet injected normally into a crossflow has 
been presented. The crossflow Reynolds number is 5600 and the blowing ratio is 0.25. Both flow 
and heat transfer predictions are presented. The following major observations are made: 

1 . Jet penetration is limited to 1 .8 hole diameters, with the primary crossflow velocity relatively 
unaffected beyong 1 hole diameter. 
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2. A weak horse-shoe vortex system consisting of two clockwise eddies is noted within one hole 
diameter upstream of the jet. This horseshoe vortex is deflected upwards and rides over the CVP 
in the near field of the jet. 

3. The near-field past jet reattachment is characterized by developing boundary layer behavior 
with ejection and sweep events commonly termed as turbulent bursts. However, these events are 
largely controlled by the large scale vortical structures. 

4. The origins of the CVP is associated with the x-vorticity generated by the interaction of the 
injected jet with the crossflow. Therefore hole-boundary-layer vorticity is not essential for the 
generation of the CVP. 

5. There is a significant amount of crossflow in the near field (x/d<5) of the jet. This crossflow 
encounters an adverse gradient along the jet-hole centerline and is responsible for the formation of 
the wall vortex system. The wall vortex is fairly stable for x/d up to 5, but beyond that is 
entrained into the CVP and is reoriented as a wake vortex system. 

6. Along the upper edges of the jet, Tollmein-Schlichting instability leads to the formation of the 
shear layer vortices. These vortices are peeled and convected downstream. 

7. The Vnns and t^ contours appear to be most strongly correlated with each other and with the 
shear layer vortices in the near field. The u^ and contours appear to be the larges tat the 
interface between the CVP and the horse-shoe vortex and along the periodic boundaries at the 
locations associated with significant spanwise crossflow. In general, > w* > v^. The 
strong correlation of these statistics with the large scale structures point to the importance of the 
latter in controlling the surface heat transfer and the near field mixing between the crossflow and 
the coolant stream. 
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8. The viscosity factors show considerably anisotropy and in the horizontal (x-z) plan show 
correlation with the large scale structures. Existing statistical models based on isotropy are clearly 
going to be inadequate in modeling this complex flow. Current efforts are focussed in examining 
the budgets of the turbulence kinetic energy and the dissipation rates, in order to provide some 
guidance toward turbulence modeling. 
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Figure 2: Comparison of DNS predictions of Pointel et al„ 1996 (— ) with the DNS predictions of 

Kim et al., 1987 ( ) and experimental data of Kreplin and Eckelman (o, +, A). Law of the 

wall shown by ( ) lines 
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Figure 3 (a) U contours at z/w = 0 (b) V contours at z/w = 0 (c) Omega X Contours at y/d 0.0239 
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Figure 4 Instantaneous velocity vectors in the cross stream planes 
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Figure 5 Instantaneous velocity vectors at (a) z/w = 0 (b) Near wall details (c) y/d = 0.0239 



Figure 6 Umean profile* at (a) r/d ■ -1 .9 (b) z/w ■ 0, Tmean profiles at (c) z/w ■ -1 .9 (d) z/w ■ 0 
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Figure 1 3 rms contours at x/d = 0.5 and 2 












Figure 1 4 rms contours at x/d = 5 and 1 0 
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Fiaure 1 5 rms contours at z/w = 
















Figure 1 7 shear stress contours at x/d = 5 and 1 0 
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Figure 18 shear stress contours at z/w 
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Figure 20 turbulent heat flux contours at x/d = 5 and 10 
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